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DYNAMICS  OF  LIQUID  FILLED  SHELL:  NON-CYLINDRICAL  CAVITY 

ABSTRACT 

A  theory  is  presented  for  the  approximate  computation  of  eigen- 
frequencies  of  liquid  oscillations  in  non-cylindrical  cavities.  The 
eigen-frequencies  are  essential  for  the  prediction  of  stability  of 
liquid-filled  shell.  The  theory  reduces  the  problem  of  finding  the 
eigen- value  to  a  simple  integration  which  can  be  performed  by  hand 
computation  when  the  shape  of  the  cavity  is  known.  Comparison  of 
theoretical  prediction  and  available  experimental,  data  shows  very 
good  agreement. 
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I.  INTRODUCTION 


The  dynamic  behavior  of  liquid-f illed,  spinning  shell  can  be  reliably- 

predicted  if  the  cavity  occupied  by  the  liquid  is  cylindrical  or  if  the 

cavity  is  spheroidal  and  completely  filled  with  liquid.  The  latter  case^ 

a  spheroidal  cavity  completely  filled  with  liquid,  is  of  relatively  little 

practical  interest  for  the  shell  designer.  In  fact  most  cavities  are 

nearly  cylindrical  and,  furthermore,  cavities  are  usually  not  100  percent 

filled.  It  is,  therefore,  very  fortunate  that  the  case  of  a  partially 

filled  cylindrical  cavity  is  accessible  to  theoretical  analysis. 

1* 

The  theory  for  this  case  was  published  by  Stewartson  in  1959 

and  has  since  become  the  most  important  tool  for  the  design  of  stable 

shell.  Although  most  practical  cavities  are  approximately  cylindrical, 

some  deviate  sufficiently  from  exact  cylindrical  shapes  that  there  is 

some  doubt  whether  treating  them  as  cylinders  is  still  a  good  approxi- 

2 

mation.  To  clarify  this  point,  Karpov  has  made  extensive  experi¬ 
mental  investigations  of  non-cylindrical  cavities.  (All  cavities 
considered  are  bodies  of  revolution. )  It  was  found  that  rounded  corners 
produce  very  little  effect  on  the  range  of  instability  but  that 
considerable  changes  result  from  modifications  like  conical  reduction  of 
one  or  both  ends  of  the  cavity.  In  view  of  the  large  variety  of  cavity 
shapes,  it  is  not  possible  to  explore  the  effect  of  cavity  shape  solely 
on  an  experimental  basis. 

* 

Superscript  numbers  denote  references  which  may  be  found  on  page  20. 
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On  the  other  hand,  an  exact  theoretical  approach  to  the  problem 
appears  to  be  quite  hopeless  for  the  following  reason:  The  -  equations 
of  the  fluid  motion  are  to  be  solved  for  boundary  conditions  that  are 
imposed  by  the  shape  of  the  cavity  walls.  If  the  cavity  is  cylindrical, 
and  only  in  this  case,  the  solution  is  separable  in  radial  and  axial  com¬ 
ponents.  The  tremendous  simplification  gained  by  the  separation  of 
variables  is  lost  if  the  cavity  shape  deviates,  however  slightly,  from  a 
perfect  cylinder.  Solutions  could  be  found,  possibly,  by  numerical., 
methods;  however,  such  an  approach  is  quite  impractical  considering  that 
the  computation  cannot  be  performed  prior  to  the  cavity  design.  The  best 
to  be  done  within  the  borders  of  an  exact  theory  is  to  compute  the  de¬ 
sired  data,  e.g.,  eigen-frequencies  -  for  a  comprehensive  class  of 
cavity  shapes.  The  shell  designer  would  then  approximate  the  actual 
cavity  by  one-  for  which  data  had  been  computed  and  tabulated.  However,  ■ 
the  advantage  of  an  exact  solution  is  lost  when  the  cavity  is  only 

approximate,  and  one  might  ask  whether  an. approximate  solution  for  the 
exact  cavity  is  not  preferable. 

II.'  THEORY  OF  APPROXIMATE  EIGEN- FREQUENCIES 
Of  particular  importance  for  the  stability  of  the  shell  are  the 
eigen -frequencies  or  frequencies  of  free  oscillation  of  the  liquid. 
According  to  Stewartson  \  instability  occurs  whenever  any  of  the 
eigen-frequencies  falls  within  a  certain  bandwidth  about  the  frequency 
of  nutation  of  the  shell.  Details  on  this  stability  theory  are  found  in 
Stewartson* s  paper.  1  Whether  a  shell  is  stable  can  be  predicted 
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when  the  eigen -frequencies  and  the  residues  of  the  forcing  term,  which 

determine  the  bandwidth,  are  known.  For  cylindrical  cavities,  the  eigen- 

frequencies  and  residues  are  computed  exactly  and  given  in  Stewartson's 
1  5 

tables.  y  For  small  deviations  from  cylindrical  shape,  small  changes 

of  the  eigen -frequencies  and  the  residues  must  be  expected.  While  accurate 

values  for  the  eigen -frequencies  are  essential  for  the  prediction  of  the 

dynamic  behavior,  the  exact  values  for  the  residues  are  less  important, 

since  these  determine  only  the  bandwidth  and  usually  the  latter  is 

affected  by  other  factors,  such  as  viscosity,  more  than  by  deviations 

2 

from  cylindrical  shapes.  From  Karpov's  experimental  data  of  un¬ 
damping  rates,  it  appears  that  it  is  sufficient  and  safe  to  assume  a 
value  for  the  residue  equal  to  that  of  an  "equivalent"  cylinder,  i.e., 
one  with  the  same  volume  and  height.  Thus,  we  are  left  with  the  determi¬ 
nation  of  the  eigen-frequencies  for  non-cylindrical  cavities.  The 
following  analysis  rests  on  the  assumption  that  the  radius  of  the 
cavity,  a,  is  a  slowly  varying  function  of  the  distance  along  the  axis, 

Z  y  1 «  6  •  ^ 


From  Stewartson's  analysis  it  follows  that  the  oscillatory  part  of  the 
pressure  -  which  results  from  the  liquid  oscillations  -  is  of  the  form: 

_/  >  i(u)t-e)  /c\ 

p  =  P(r,z )  e  (2) 

where  P(r,z)  satisfies  the  differential  equation: 
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(3) 


with 


a2p  i  ap  i 

3ra  r  clr  r3 


P  =  Q3 


a2p 

5F 


-  k 
(l  ~-~t)2 


1 


(4) 


tu  is  the  eigen-frequency,  (r,  9,  z)  are  polar  coordinates  and  t  is 

ID 

the  dimensionless  eigen-frequency,  t  *  —  ,  where  fi  is  the  frequency  of 
spin . 

Let  us  use  the  notation  (u,  v,  w)  for  the  components  of  the  oscillatory 

velocity  in  the  polar'  coordinate  system  (r,A,z).  Then  (u,  v,  w)  can  he 

expressed  by  certain  linear  combinations  of  the  pressure  and  its  partial 

derivatives  and 

°r  3z  , 

We  assume  that  the  cavity  shape  is  given  by: 


a  =  a(z)  ,  0  s  2  s  2c 


(5) 


and  that  Equation  (l)  holds  for  0  ^  z  <  2c  .  If  the  cavity  is  partially 
filled  with  liquid  it  is  assumed,  as  in  Stewartson's  case,  that  the' 
undisturbed  free  surface  is  a  cylinder  of  radius  b  .  For  the  following 
analysis  we  must  assume  that 

b  <  a(z)  for  0  s  2  s  2c  (6) 

Equation  (3)  must  be  solved  in  connection  with  certain  boundary  conditions, 

e  S  •  5 

u  =  t—  w  at  r  =  a(z)  (7) 


10 


w  =  0  at  z  =  0,  2c 


•  (8) 


and  one  boundary  condition  at  r  =  b,  which  requires  that  the  pressure  be 

zero  on  the  free  surface  when  the  cavity  is  partially  filled.  For  a 

better  understanding,  a  brief  review  of  the  eigen-value  problem  for 

Stewartson's  case  is  given  in  the  following:  For  the  cylindrical  cavity, 
da 

i.e.,  if  —  =  0,  the  solution  of  (3)  with  appropriate  boundary  con¬ 
ditions  can  be  found  by  separation  of  variables: 

P  =  C(okr)  .  cos  kz  (9) 

where  C  =  A*  Jj.  (okr)  +  B’Yi(okr)  and  Of  ,  Yj  are  Bessel  functions  of  the 
first  and  second  kind  with  the  argument  okr  .  The  solution  (9)  must 
satisfy  certain  boundary  conditions  at  r  *  b  and  r  =  a  .  Without 
giving  a  detailed  derivation  (which  can  be  found  in  Reference  l)  it  suffices 
here  to  state  the  following  results:  For  given  ba/as  and  t  the  boundary 
conditions  at  r  =  b  and  r  =  a  lead  to  a  transcendental  equation  for  ka 
with  discrete  roots: 

ka  =  (b2/a3  ,  t)  ,  n  =  1,  2,  3,  ....  (10 ) 

where  n  is  related  to  the  number  of  radial  waves  of  the  solution  (9). 

The  boundary  conditions  at  z  =  0  ,  2c  lead  to  a  similar  equation 
for  kc  with  roots: 

kc  =  |  [2J+1]  0=0,  1,  2  ...  (11) 
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Since  w  is  proportional  to  ^  ,  Equation  (8)  can  be  written: 

z  -  0,  2c  .  (12) 


kc  =  |  .  m  •  (13) 

Because  of  certain  symmetry  requirements,  the  number  of  half-waves  in  the 
z-direction  must  be  odd,  i.e.,  m  =  (2j+l)  and  Equation  (ll)  follows  from  (13). 
Comparison  of  (ll)  and  (10)  gives  a  condition  for  the  fineness  ratio: 

c  tt/2  (l4) 

a[2j+l]  =  \  (b2/aa  ,r) 

Q 

Equation  (l4)  expresses  that  only  certain  discrete  fineness  ratios  — 

exist  for  given-  bs/a3  and  t  or,  for  given  b3/a3  and  c/a  there  exists 

a  set  of  discrete  eigen-functions  TJn  according  to  the  choice  of  j  and  n 

1  c 

in  Equation  (l4).  In  Stewartson’s  tables  the  various  a(2j+i)  are 
tabulated  as  functions  of  b2/a3  and  t  .  Let  us  call  these  functions 
S„  (b2/a3  ,  t)  ,  i.e. , 

fej+l)  =  ^  (b2/a2  ’  t)  =  T^(b3/a3  ,  t J 


Equations  (12)  and  (9)  yield: 


sin  .k  2c  =  0  , 


(15) 


Returning  to  the  eigen-value  problem  for  the  non-cylindrical  cavity,  we 

notice  that  the  only  difference  is  in  the  boundary  condition  (7).  This 

condition  states  that  the  velocity  component  normal  to  the  wall  r  =  a 

is  zero.  If  now  the  inclination  of  the  wall  toward  the  axis  is  small,  i.e., 

|~|  4  1  ,  an  obvious  approximation  of  Equation  (7)  is: 
u.z 

u  =.  0  at  r  =  a(z)  (16) 

The  form  of  the  approximate  boundary  condition  (16)  suggests  searching  for 
approximate  solutions  of  Equation  (3)  which  are  of  the  same  general  form  as 
Equation  (9)  except  that  the  radial  wave  number,  k,  should  be  a  slowly 
varying  function  of  z  according  to  the  fact  that  the  radius,  a,  changes 
slowly  with  z.  Thus,  we  try  as  an  approximate  solution  of  (3): 

P  =  C(o'kr)  .  cos  $(z)  (17) 

where  C  is  a  cylinder  function  of  the  argument  ukr  and  k  depends 
weakly  on  z  .  The  phase  $(z)  can  no  longer  be  assumed  to  be  equal  to 

rz 

kz  but  rather  $  =  kd£ 

o 

It  can  be  verified  easily  that  (17)  approximately  satisfies  the 
differential  Equation  (3)  provided  that: 

(18) 

and 

IH1  <  &  (19) 

~  '  t  •  "1 

AEHaiUu &  1.  .  .  ;■  UD< 

. . 
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.  Since  the  boundary  condition  at  r  =  b  and  the  approximate  boundary 
condition  at  r  =  a  (Equation  (l6))  are  -  locally  -  the  same  as  for  the 
cylindrical  cavity,  we  arrive  naturally  at  the  same  condition  for  ka  , 
viz.,  according  to  Equation  (10), 

k  =  7  \  (Wa2,  t)  (20)- 

Equation  (20 )  gives  the  z-  dependent  wave  number  k(z)  and,  at  the  same 
time,  shows  that  the  condition  .  (19)  is  a  consequence  of  |^|  4  1  (since 
and  its  derivative  with  respect  to  ba/as  are  of  order  unity). 

Thus,  we  conclude  that  (17)  with  (18)  approximately  satisfies  the 
differential  Equation  (3)  and  the  boundary  conditions  at  r  =  b  and  r  =  a. 

Finally,  the  boundary  condition  (12);  =  q  at  z  =  0,  2c 

requires  that  sin  $  =  0  at  z  =  0,  2c  or: 

§(o)  =  0  (21) 

2c  ’ 

$(2c)  =  f  k(z)  dz  =  n[2j+l]  (22) 

0 

Equation  (22)  can  be  considered  as  a  generalization  of  Equation  (11). 

Substituting  (20)  into  (22)  gives: 

2c 

tt[2j+1]  =  f  ~  \(^8A2,  t)  dz  (23) 

cl 

0 

Equation  (23.)  connects  b  and  t,  i.e.,  if  the  radius  of  the  cylindrical 
void,  b,  is  given,  (25)  determines  t  and  vice  versa. 


According  to  Equation  (15)  the  function  T|n  is  related  to  Sn  ,  a  function 
which  is  tabulated  in  Stewartson's  tables.  Thus,  it  is  convenient  to 
rewrite  (23)  by  substituting  according  to  Equation  (15)  into  (23)  and, 
after  dividing  by  rr(2j+l),  one  obtains: 


2c 


1  = 


1 

2c 


sMZi+ll  ... 

S,,  (b2/a2 ,  t) 


dz 


(24) 


Equation  (24)  is  easy  to  memorize:  Sn(bs/a2,  t)  is  just  the  c/a(2j+l) -value 
of  a  cylinder  with  eigen-frequency  t  and  fill-ratio  b2/as  .  Thus,  the 
integrand  of  (24)  is  the  ratio  of  the  local  c/a(2j+l)  -  value  and  the 
c/a(2j+l)  -  value  which  would  correspond  to  t  and  the  local  b2/a3  -  value. 
Equation  (24)  then  states  that  the  mean  value  of  this  ratio  -  averaged 
over  the  length  of  the  cavity  -  should  be  equal  to  one . 

Usually,  the  higher  radial  modes  are  unimportant,  so  that  (24)  must 
be  solved  only  for  n  =  1  .  We  will  therefore  write  just  S  instead  of 
Sn  ,  having  in  mind  that  (24)  is  valid  for  any  of  the  S„  values. 

The  evaluation  of  Equation  (24)  is  difficult  when  S  is  given 
numerically  (as  in  Stewartson's  tables),  since  one  has  to  assume  both, 
b  and  t,  in  order  to  perform  the  integration  numerically  and  eventually 
repeat  the  procedure  with  changed  values  of  t  (or  b)  until  the  correct 
value  of  t  (or  b)  can  be  observed  by  interpolation.  For  a  100  percent 
filled  cavity,  i.e.,  b  =  Q,  the  integration  simplifies  considerably, 
since  S  becomes  independent  of  z  and  can  be  taken  out  of  the  integral. 

One  obtains  from  (24): 
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2c 

.  s<0’  T)  ■  2c  I  afsjrry  iz  (25) 

0 

C 

Since  S  is  the  a(gj+i)  “  value  of  a  cylinder  with  eigen-frequency  t, 
Equation  (25)  can  be  interpreted  in  the  following  way:  A  completely  filled 
non-cylindrical  cavity  has  the  same  eigen-frequencies  Tn  s  as  an  "equivalent 
cylindrical  cavity",  which  is  defined  by  having  a  fineness  ratio,  c/a  equal  to 
the  averaged  c/a  of  the  non-cylindrical  cavity.  For  a  100  percent  filled 
cavity,  therefore,  only  one  integration  must  be  performed  to  determine  the 
average  c/a. 

One  way  of  solving  the  eigen-value  problem  given  in  Equation  (24)  is 
to  approximate  l/S  by  a  power  series  in  b2/a2  .  If  T  is  given  and  b 
is  to  be  determined  one  could  plot  l/S  ,  for  the  particular  t  given, 
versus  b2/a2  and  approximate  the  obtained  curve  by  a  polynomial  in  b2/a2 , 

If  neither  T.  nor  b2/a2  are  too  large,  the  following  formula  is  convenient, 
and  quite  accurate: 

I  ‘  i  Ti +  ^  (?)’]  ’  <26> 

S  is  the  value  S(o,  t)  as  obtained  from  Stewartson's  tables  for 
o 

b2/a2  _=  0,  The  approximation  (26)  is  valid  within  the  following  limits: 

0  S  t  £  0.12 

'  0  s  b3/a2  £  0.15  .  ^ 

with  (26)  substituted  into  (24)  one  obtains: 
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(28) 


2c 

so(T)  =  2^  J  a(2J+l) 
0 


[ 


1  +  1.2' 


dz 


Equation  (28)  has  the  advantage  that  the  right  side  is  independent  of  t 
and  the  left  side  independent  of  b,  i.e.,  by  integration  of  (28)  one  finds 
a  relation  of  the  form: 


S„(t)  =  Cj  +  c2  b4  (29) 

A  better  approximation  than  (28)  may  be  obtained  by  choosing  a  coefficient 
different  from  1.26  in  Equation  (26)  adapted  to  the  particular  t,  if  t  is 
known.  One  could  also  include  higher  powers  of  b2/a2. 
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III.  COMPARISON  WITH  EXPERIMENTS 

In  order. to  test  the  theory  given  in  the  preoeeding  section  some  of 

'  p 

the  experimental  findings  of  Karpov  have  been  compared  with  theoretical 
predictions.  The  cavities  investigated  by  Karpov  were  composed  of  a 
cylindrical  and  a  conical  section.  A  schematic  of  a  cavity  shape  is  shown  in 
the  following  figure.  In  the  diagram,  the  filling  ratio  which  corresponds 
to  a  fixed  eigen-frequency  t0  is  plotted  versus  the  variable  cone-angle  fl. 
Three  curves  are  shown  that  correspond  to  three  different  ratios. of  conical 
section  h  to  overall  length  2c  . .  For  9=0  all  cavities  degenerate  into 
a  cylinder  of  fineness  ratio  —  =  2.687  and  the  filling  ratio  (l  -  bs/a2 ) 

SL 

attains  Stewartson's  value,  (it  is  j  “  1,  n  =  l).  The  solid  curves  in  the 
figure  show  the  theoretically  predicted  values,  the  symbols  the  experimental 
data.  The  theory  obviously  gives  good  predictions  even  for  0 -values  for 
which  the  supposition  on  which  the  theory  rests,  viz.  1^1  =  4  1  ,  is 
hardly  fulfilled.  The  reason  for  this  is,  roughly,  that  the  neglected  terms 
in  the  boundary  conditions  and  in  the  differential  equation  are  orthogonal, 
or  nearly  so,'  to  the  other  terms.  Therefore,  although  the  "local  error"  is 
of  order  j~^|  ,  the  "averaged  error"  over  the  length  of  the  cavity  can  be 

»  CLZ 

very  small.  The  theoretical  curves  in  the  figure  were  obtained  by  solving 
Equation  (25)  for  fixed  t  ~  0,055-  To  perform  the  integration,  the 
following  approximation  was  used: 

?  -  0.944  +  1.074  (g  )* 

The  integration  could  then  be  performed  analytically  .and  b2/a2  is 

*  h.  * 

found  as  a  function  of  0  and  rr~ 

.  2c 
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Filling  Ratio  At  Fixed 
Resonance  Frequency  T0=Tn  =  0.055 


9  -  Degrees 
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